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ABSTRACT 

Observations of prominences show them to require well-developed magnetic shear 
and to have complex small-scale structure. We show here that these features are 
reflected in the results of the theory of radiative condensation. We have studied, 
in particular, the influence of the nominally negligible contributions of 
perpendicular (to B) thermal conduction. We find a large number of unstable modes, 
with closely spaced growth rates. Their scale widths across B show a wide range of 
longitudinal and transverse sizes, ranging from much larger than to much smaller 
than the magnetic shear scale, the latter characterization applying particularly in 
the direction of shear variation. 


INTRODUCTION 

Coronal prominences owe their existence to a condensation process which occurs 
as a consequence of an instability in the thermal equilibrium of a diffuse medium 
(Parker, 1953; Field, 1965; Hildner, 1974). The condensation mechanism relies on 
optically thin radiation whose dependence on thermodynamic variables (e.g., density 
and temperature) is such that a cool, dense perturbation loses more energy through 
radiation than it gains through adiabatic and non-adiabatic heating processes and 
thermal conduction. 

Prominences and filaments (as seen on the disk) in the solar atmosphere are 
often observed to form above a magnetic neutral (polarity-inversion) line in regions 
of Increasing magnetic field shear as indicated by photospheric magnetograms 
(Martin, 1973; Leroy, 1978). Local heat conduction, which is strongly attenuated in 
directions perpendicular to the magnetic field, would dominate radiation and other 
forms of energy transport if no field were present, thereby suppressing the thermal 
instability. Thus, one can expect the equilibrium field structure to exert a strong 
influence over the formation of prominences. 

These empirical and physical considerations have motivated previous 
computational studies of the dynamics of the thermal instability in a sheared 
magnetic field (Chiuderi and Van Hoven, 1979; Van Hoven and Mok, 1984; Van Hoven et 
al., 1984; Sparks and Van Hoven, 1985). One result of these theoretical studies is 
that a sheared background field is necessary for the existence of a true localized 
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condensation* Secondly, it is surprising that the width of the condensation in a 
direction perpendicular to the shear layer does not correspond to those points at 
which the magnetic field has tilted sufficiently that the radiation loss is roughly 
balanced by parallel thermal conduction, nor is it affected by realistic values of 
the perpendicular conductivity. 

The present study is devoted to the delineation of some additional consequences 
of the presence of anisotropic thermal conduction in a sheared-field filament* One 
finds that new unstable excitations appear, with complicated transverse variations 
(multiple nodes), and that perpendicular conduction provides modes with shorter 
wavelengths and faster growth. 

FORMULATION 

To model a sheared, active-region, magnetic field, we use the planar 
form B/B = F(y/a)e + G(y/a)e where F(0) « 0 and F^ + * 1, which is 

with uniform temperature (T) and density (p). [The example we use is 
F(y/a) ■ tanh y/a.] 

We describe the plasma dynamics by the compressible ideal MHD limit, 
energetics by the heat-transport equation (Chiuderi and Van Hoven, 1979) 

_ ypV»v - (y-l)[H o - p 2 $(X) + V«jc*VT] (1) 

which includes an unspecified (and unknown) constant heat input for thermal- 
equilibrium balance, optically thin radiation losses and anisotropic heat flow. If 
one considers T(y,z,t) = Tq + T^(y)exp(vt + ikz), one can linearize (1) as 

5T dp 

-T-— + (y-1) — — = Q T - Q n, 

dt w ' dt pi T H 1 

- [k 2 a 2 F 2 Q (| + k 2 a 2 G?Q - Q (Tj/Tj)] (2) 


force-free 

consistent 


and the 


where T^ and p, are fractional perturbations and TjT * adT^/dy. The heat-flow rates 
for the non-adiabatic terms on the right side of (2) include the radiation rate 
Q at constant density and^the generalized parallel-plus-perpendicular thermal 
conduction rate Q « icT^/a p Q in the square bracket on the right. 

When the two-dimensional dynamic equations are simplified, they reduce to the 
set of coupled equations (Van Hoven and Mok, 1984; Sparks and Van Hoven, 1985) 

.. 2 „ 2 ' 2 2 W 2 . 2_,2 " 1 ., 2 . 2 . 
q - a F (v +a )(v +a F ) q - (a +a^ )q = 0 

-1 9 9 <3) 

T x = [(Y-Dv - 2 T ][v-Q p -K J k ] 1 («yv )q 

where q = is the total pressure perturbation. The normalized wave 

numbers are a = ka anS 
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a 2 = v 4 (v - Q + Q ) 

V p K 

X [ V 2 ( v-Q p -Ky + V 2 PY ( v 2 +a 2 F 2 ) ( v-Qp+Q^/y ) ] _1 


which depends implicitly on T^"/T^ (through Q^) when 0 and provides the only 
energy-transport contribution to the _q equation. 

The solutions of equations (3), with boundary conditions requiring localization 
of the excitation near y = 0 and exponential decay as y -► ± «, provide 
eigenf requencies (growth rates) v which depend on ka, and eigenfunctions T^(y) # 
Significant information about the allowable values of v(ka), and about the 
structure of the eigenfunctions, can be obtained from the poles and zeros of a and 
from considering local solutions in a uniform field (Chiuderi and Van Hoven, 1^79; 
Sparks and Van Hoven, 1985). 

RESULTS 

The simplest radiative-instability case to consider in a sheared field has 
Q = 0 so that radiation merely competes against adiabatic compression. The 
solutions of (3a) then provide a series of modes which have the essential 
characteristics of the solutions to more complete formulations. The modes exhibit 
increasing numbers of y-direction nodes, and growth rates approaching 
Q = (Q + V 2 PyQ )/(!+ V 2 py), the radiation rate for which perpendicular (to B) 
piasma Motions oEcur at constant total pressure. In fact, these solutions have the 
typical property that the total-pressure perturbation is much smaller than the 
thermal pressure p^, especially at shorter wavelengths (Sparks and Van Hoven, 1985). 

The addition of the nominally dominant effects of parallel thermal conduction 
does not change the situation very much. The principal modification is the 
introduction of a pole in a ,^which prevents the existence of solutions at 
wavelengths shorter than the s 4 * ® curve of Fig. 1 (Chiuderi and Van Hoven, 

1979). 

It is necessary, finally, to add the effects of perpendicular thermal 
conduction Q , which should not (a priori) be important, to be able to 
obtain a reasonably complete treatment of the structure and growth of these 
sheared-f ield radiative modes (Van Hoven and Mok, 1984; Van Hoven et al., 1986). 

In order to see the important effects of Q^, I have shown a qualitative growth- 
rate curve in Fig. 1, which also displays the various instability rates identified 
by Field (1965). The behavior of the eigensolutions is different on the two sides 
of the curve s Z = a +a (F-1,Q *0) -► » [or a ■ a ( v ) ] . The growing modes on the 
lower left are variants of the original sheared-?ield modes found by Chiuderi and 
Van Hoven (1979). The addition of k to the energy transport resolves the steep 
gradients of these solutions (Van Hoven and Mok, 1984) and allows them to have 
multiple radial nodes, within a width given by Eq. (18) of this earlier paper. It 
Is somewhat unusual that the "fundamental" (no nodes) transverse-variation mode has 
the lowest growth rate. The practical consequence of this fact is unclear, however, 
since the growth-rate curves effectively lie on top of each other for typical solar 
coronal parameters • 
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Figure !• Radiative f ilamentation growth rates vs wavenumber. 


The faster-growing, shorter-wavelength modes on the upper right in Fig. 1 only 
appear as a result of the presence of k . (Van Hoven, Sparks, and Tachi, 1986). 

[Small as this coefficient is (k /k^ ~ 10 in the corona), it resolves certain 

singularities which appear when £ 0.] For these modes, the fundamental is the 
fastest growing, although the inter-mode spacing near the peak is (again) negligible 
in practice. The transverse structure of these eigensolutions is unusual, for 
v > Q , in that the (negative) peak of the temperature perturbation is located away 
from ?he axis on each side (Van Hoven, Sparks, and Tachi, 1986). The inner and 
outer edges of the temperature peak(s) can be scaled from the equivalent of Fig. 1. 
One takes the a positions [a (v)] of the s^ = <» curve at the relevant eigenvalue of 
v(a), divides it by a - ka afid calculates y/a - tanh -1 (ct /a)» The resulting widths 
for reasonable parameters are ~ 10 a (located at ~ ± P 10~a) for X/2 ~ 0.3a, but 
one must remember that a number of modes grow at nearly the same rate. The outer 
edges of these modes, for which a ® y(Q -v)/Q , exhibit the only direct dependence 
on the parallel conductivity. p p 
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DISCUSSION 


The principal conclusion of our study of the linear eigensolutions of the 
radiative filament-condensation instability in a sheared magnetic field is that the 
theoretical results are nearly as complex as the observational results. The range 
of 


growth ratesfor the radiative instability is small, with - 

U ~ 10 - ^ * n. T, Z < v < Q - 1.8Q for the usual estimate $(T) « T -1 (Hildner, 1974) 
p 9 6 ~ p p 

6 9-3 

at coronal temperatures T ~ 10 T^ K and number densities n ~ 10 n^ cm . 


Not only is the range of growth rates narrow, but there are a large number of 
distinct eigenmodes within this range. There are two groups, one with wavelengths 
(the vertical width of the characteristic knife-blade form) greater than the 
magnetic shear scale, and the second with mainly shorter wavelengths. In most 
cases, there is an approximate balance between magnetic pressure increases and 
thermal pressure decreases. 


These eigenmodes also exhibit a complex structure in the transverse direction, 
equivalent to the horizontal thickness of the knife-blade filament. The temperature 
(and density) profiles oscillate in this direction, with a number of nodes. The 
long wavelength modes are concentrated in the center of the shear layer, but the 
shorter modes often show a hollow profile with the coolest layers separated from the 
shear center. These latter excitations extend to the point where radiation is 
overtaken by parallel thermal conduction (which is relatively strong for short 
wavelengths) . 


We will not know which of these many excitations is (are) the dominant one(s) 
until we complete a series of nonlinear computations which are now in progress. 


I wish to acknowledge the many contributions of Drs. L. Sparks and T. Tachi to 
this investigation. This work was supported, in part, by the National Science 
Foundation and by the National Aeronautics and Space Administration. 


REFERENCES 


Chiuderi, C. and Van Hoven, G. : 1979, Astrophys. J. Letts. 232 , L69. 

Field, G. B.: 1965, Astrophys. J. 142, 531. 

Hildner, E.: 1974, Solar Phys. 35, 123. 

Leroy, J.-L.: 1978, Astron. Astrophys. 247. 

Martin, S. F.: 1973, Solar Phys. 31 , 3. 

Parker, E. N. : 1953, Astrophys. J. 117, 431. 

Sparks, L. and Van Hoven, G. : 1985, Solar Phys. 97 , 283. 

Van Hoven, G. and Mok, Y.: 1984, Astrophys. J. 282, 267. 

Van Hoven, G. , Tachi, T., and Steinolfson, R. S.: 1984, Astrophys. J. 280 , 391. 
Van Hoven, G., Sparks, L. , and Tachi, T.: 1986, Astrophys. J. 300, 249. 


